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^ : Abstract 

,^ ' The strong phases and CP violation in the rare B^ — > K~^K~, K K^, K*^K*~ decays are inves- 

ji. . tigated. As these decays proceed only via annihilation type diagrams in the standard model (SM), 

(— I , a dynamical gluon mass is introduced to avoid the infrared divergence in the soft endpoint regions. 



The Cutkosky rule is adopted to deal with a physical-region singularity of the on mass-shell quark 
propagators, which leads to a big imaginary part and hence a large strong phase. As a consequence, 
large CP asymmetries are predicted in those decay modes due to a big interference between the an- 
nihilation amplitudes from penguin and tree operators, which may be tested in future more precise 
experiments. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.38.Bx 



I. INTRODUCTION 

Charmless i?-meson decays are of crucial importance to deepen our insights into the flavor 
structure of the standard model (SM), the origin of CP violation, the dynamics of hadronic 
decays, and to search for any signals of new physics beyond the SM. The CP violation, as an 
important part of B physics, has been paid much attention in recent years. Both mixing-induced 
and direct CP violations have been observed in the neutral B meson decays |]J, 12, S 1^ , |5| , which 
has provided a new window for exploring the origin and mechanism of CP violation after the 
estabhshment of indirect and direct CP violations in kaon decays la |7|, |8|, |9|, |lO| . The possible 
implications of charm! 



recent papers 
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ess B decays and their large CP violation have been investigated in the 
131 . In the SM, the only source of CP violation in the SM is the single 
Kobayashi-Maskawa phase |lj| in the mixing matrix that describes the charged current weak 
interaction of quarks. However, physics beyond the SM are usually associated with new sources 
of CP violation. For instance, rich sources of CP violation can be induced from a single relative 
phase of vacuum in the simple two Higgs doublet model with spontaneous CP violation(S2HDM 
or type III 2HDM) |l^. The model can provide a natural explanation for the CP violation in 
the SM and also lead to a new type of CP- violating source with each quark and lepton carrying 
a nontrivial CP-violating phase. Therefore, explorations of CP violation may well indicate 
physics beyond the SM, or may be very helpful to distinguish between various realizations of 
one particular kind of new physics after the corresponding new physics particles have been 
observed directly. 

In addition to being served as a tool of looking for any new physics, studying of CP vi- 
olation can also be used to test various factorization hypotheses, such as the "naive" factor- 
ization approach (FA) [la], the QCD factorization (QCDF) (l2|, and the perturbation QCD 
method (pQCD) [18[. These methods have very different understandings of hadronic i?-meson 
decays: For the FA method, it cannot predict the direct CP asymmetries properly due to the 
assumption of no strong re-scattering in the final states; the pQCD generally predicts large 
strong phases and big direct CP violations. Furthermore, this approach also believes that an- 
nihilation diagrams are important in non-leptonic two-body 5-meson decays (l9l.l2fll]: while the 
QCDF favors small direct CP violations in general because of the a^-suppressed strong phases. 

It is known that, in most cases of two-body 5-meson decays, the weak annihilation contri- 
bution carries different weak and strong phases for the tree and penguin amplitudes, which is 



very important for studying CP-violating observables. Meanwhile, the calculation of annihi- 
lation contributions is interesting by itself, since it can help us to understand the low energy 
dynamics of QCD involved in heavy meson decays and the viability of the theoretical ap- 
proaches mentioned above. Motivated by these arguments, we shall investigate in this paper 
the pure annihilation processes B^ -^ K^K~, K*^K^, K*^K*~ . These decay channels have 
some interesting features: firstly, they are all pure annihilation processes and studying these 
decay modes in the SM can serve as a probe for the annihilation strength in hadronic 5-meson 



decays; secondly, the current experimental data on the decay B^ -^ K^K~ 2l| has already 
provided a direct evidence for the existence of the annihilation contributions. By comparing 
our theoretical result for the branching ratio of the decay B^ -^ K^K~ with the experimental 
data, we can test the feasibility of the theoretical method and get a deeper insight into the 



penguin annihilation and the VT-exchanged topologies in B ^ txtx ^iiK decay modes [13|, |22| |: 
thirdly, their strong phases are calculable and found to be big, which leads to large CP violation. 
The resulting branching ratios for B{B^ -^ K'^K") is consistent with the current experimental 
data [2]J; Finally, studying the polarization in the 5° — *> K*^K*~ decay may also help us to 
clarify whether annihilation contribution could resolve or alleviate the polarization puzzles in 
B — > (pK* decays as suggested in Ref. 



In calculating the hard scattering kernel, we shall use the Cornwall 2J| prescription of gluon 
propagator by introducing a dynamical mass of gluon to avoid enhancements in the soft endpoint 
region. It is interesting to note that recent theoretical [25| and phenomenological J26| studies are 
now accumulating supports for a softer infrared behavior for the gluon propagator. Moreover, 
we will adopt the Cutkosky rule J22| to deal with the physical-region singularity caused by the 
on mass-shell quark propagators, which then produce big imaginary parts from the kinematic 
region where the virtual quark becomes on mass-shell. By applying these two "tricks", we 
observe that the main contributions to the decay amplitudes come from the nonfactorizable di- 
agrams, and the CP-averaged branching ratios of 5° -^ K^K~ ,K^K*~ , K~K*~^, and K*'^K*~ 
are estimated by using the QCD factorization to be about 2.02 x 10"^, 4.23 x 10"^, 5.70 x 10"^, 
and 6.89 x 10^*^ for a given gluon mass rrig = 500 MeV, respectively. Moreover, big strong 
phases are predicted in these decay modes, and hence the direct and mixing-induced CP vio- 
lations Ckk and Skk are found to be very large in these decays as the differences ACkk and 
ASkk are small in B^ -^ K^K*^ decays. The predictions may be tested in the more precise 



experiments at B-factory and the LHC-b. 

The paper is organized as follows. In section II, we first analysis the relevant Feynman 
diagrams and then outline the necessary ingredients for evaluating the CP asymmetries of the 
B^ -^ K^K~ , K*^K^, K*^ K*~ decays. In section III, we present the approaches for dealing 
with the physical-region singularities of gluon and quark propagators. The numerical results 
of the CP-averaged branching ratios and large strong phases are given in section IV. Finally, 
we discuss CP asymmetries for those decay modes and give conclusions in sections V and VI, 
respectively. The necessary input parameters and the complete decay amplitudes for those 
decay modes are given in the appendixes. 

II. REPHASE-INVARIANT CP- VIOLATING OBSERVABLES 

Using the operator product expansion and renormalization group equation, the low energy 
effective Hamiltonian for charmless two-body i?-meson decays can be written as J2q] 

where Ci{fi){i = 1, ■ ■ ■ , 10) are the Wilson coefficient functions which have been reliably eval- 
uated to next-to-leading logarithmic order. The effective operators Oi are defined as follows: 

Oi = {diUi)v-.A{Ujbj)v-A, O2 = {diUj)v-A{Ujbi)v-A, 

O3 = idibi)v-A Y.iQjQj)v-A, O4 = {dibj)v-A Y.i9j9i)v-A, 
q q 

O5 = {dibi)v-Aj2iQjQj)v+A, Oq = {dibj)Y_Aj2i9j'li)v+A, C2) 

q q ^ ^ 

O7 = ^{dibi)v-AY.'^qiQjQj)v+A, Os = ^{dibj)v-AY. eq{qjqi)v+A, 
Q q 

O9 = ^idibi)v-AY.'^qiQjQj)v-A, Oio = l{dibj)v-AY.eq{qjqi)v-A- 
q q 

Here i and j are SU(3) color indices, q denotes all the active quarks at the scale fi = 0{mb), 
i.e., q = u,d, s, c, b. 

With the effective Hamiltonian, calculations of the leading order amplitudes for B^ — > 
K~^K~ ,K*^K'^,K*"^K*" decays are straightforward. However, due to the conservation of 
the vector current and partial conservation of the axial- vector current, the leading order am- 
plitudes will vanish in the limit m^, m^ — *> 0. In order to probe the annihilation strength and 



discuss CP violation in these processes, we have to consider the next-to-leading order {as order) 
contributions. 

Up to the as order, the relevant Feynman diagrams contributing to the B^ -^ K^K~ 
decay (the corresponding diagrams for K*^K^ and K*^K*~ modes are the same) are depicted 
in Fig. n (corresponding to the ly-exchanged annihilation diagrams) and Fig. |21 (corresponding 
to the penguin annihilation diagrams and including the case with the exchange of m v-> s). For 
the factorizable diagrams (a) and (b) in Figs. ^ and |21 their contributions cancel each other 
both in the PF-exchanged and in the penguin annihilation diagrams, so that the non-factorizable 
contributions will dominate the decay, which can be obtained by calculating the amplitudes of 
diagrams (c) and (d) in Figs. Q and El Moreover, it is noted that these properties hold in all of 
these decay modes. 

The direct CP violation occurs only if there are two contributing amplitudes with non-zero 
relative weak and strong phases. The weak phase difference can arise from the interference 
of amplitudes from various tree(current-current) and penguin diagrams. From the Feynman 
diagrams in Figs.Hand|2l we can see that in these decays there are two kinds of CKM elements, 
VubV*^ from tree operators and VtbV^*^ from penguin ones, which will induce weak phase difference 
and hence CP violation. The total decay amplitudes for B^{B^) — > K^K~ mode can be written 
as 

^(50 ^ j^.+^.-) ^ V:,VudAT - VlVuAp = V:,V^,At[1 + ze'^''+% (3) 

A{B' ^ K^K~) = VutV:,AT - VtbV;,Ap = VutV:,AT[l + ze'^--+% (4) 

where a = arg[-V;^\4d/\/*^Kd], z = \VtlVtd/V*^Vud\\Ap/AT\, which indicates the interfer- 
ence strength between the annihilation amplitudes from penguin and tree operators, and 
S = arg(^p/^7^) is the relative strong phase between the penguin and the tree annihilation 
amplitudes. A similar consideration can be applied to the B^ -^ K^K*^ and B^ -^ K*^K*~ 
decays. The resulting decay amplitudes by using QCD factorization are given in Appendix B. 
For neutral 5-meson decays, the time-dependent CP asymmetries are defined as 

A ... _ r(i?°(t) ^ /) - T{B\t) -. /) 

^cp\ ) Y{B%t) ^ /) + r(50(t) ^ /) ■ ^ ' 

When the final state is a CP eigenstate, such as the B^{B^) — > K^K~ decay, the time- 
dependent CP asymmetries can be written as 

Acpif) = CKKCOs{Amt) + S'/<xsin(Amt), (6) 







FIG. 1: The W^-exchanged annihilation diagrams for B^ — > K~^K decay. 







FIG. 2: The penguin annihilation diagrams for B^ — > K'^K decay. 



where Am is the mass difference of the two eigenstates of Bd mesons. Ckk and Skk are 
parameters describing the direct and the mixing- induced CP violation, respectively. Both of 
them depend on the CKM and hadronic matrix elements 

— 2Jm(Acp) 



C 



KK 



1 I \ |2 

1 + |AcpP' 



S 



KK 



l + |Ac"i2 



\CP| 



with 



A 



CP 



V,lVtd{K^K~\H,jf\B^) 



(7) 



(8) 



From Eqs. Q and (JH), the CP- violating parameters Ckk and Skk can be expressed explicitly 



as 

— 22; sin a sin(5 — sin 2q; — 22; sin a cos 5 . , 

^KK ^ 1 , o T"^ — ^) '-'^i^' = "TT^ x~, — 2~' ^^) 

1 + 2z cos a cos + z'' 1 + 2z cos a cos + z'^ 

which shows that both the direct and the mixing-induced CP violation depend not only on 
the strong phase 6 but also on the magnitudes of z. Thus, when the contributions of different 
weak decay amplitudes are comparable each other, there will be a high likelihood for observable 
CP-violating asymmetries. 

For the case of -B ^ K*^K^ decays, as the final state is not an CP eigenstate, the CP- 
violating asymmetries become complicated. There are in general four decay amplitudes which 
can be expressed as 

g = {K+K*-\H,ff\B), h = {K+K*~\H,ff\B), 

g = {K-K*+\H,ff\B), h = {K-K*+\H,ff\B). 



Following the discussions in Ref. |22(, there exist in general four rephase-invariant parameters 
a^/, a^/, a^+e', fle+e' which characterize the CP asymmetries. We may redefine the following four 
parameters 



CkK — 2{(^e' + Of) 5 ^CkK — 2{0't' " 0,1' 

SkK = nide+e' + Cte+e')' ^SkK = 9 («£+£' — C 



(11) 



with 



\9?-\h? _\W-\9? 

Ojfi — -; — 77. 1 . ,„ , dfi 



-2Im(V^) -2lm{g/h) 



(12) 



'" 1 + iv^r " i + \9ih\^ 

For final state being CP eigenstate, one has ACkk = and ASkk = 0. 

As for the B -^ K*^K*~ decay mode, since the total amplitudes are dominated by the 
longitudinal ones, which can be seen below, one can evaluate its CP asymmetries in the same 
way as for the B —^ K^K~ decay. 

III. TREATMENTS FOR PHYSICAL-REGION SINGULARITIES OF GLUON AND 
QUARK PROPAGATORS 

To perform a numerical calculation, the QCD factorization approach may be used to eval- 
uate the amplitudes of B^ -^ K^ K~ , K*^ K^ , K*^ K*~ decays. The details are presented 



in Appendix B. In Eqs. ()24|) - (f^ . one will encounter the endpoint divergence, which is the 
most difficult part to deal with in the annihilation diagrams within the QCD factorization 

Xa =_L1 ±_PAe'^) In ^ in the lit- 



framework. Instead of the widely used treatment J^ -^ 



erature 
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32], we shall use an effective gluon propagator J2J, |33| to treat the infrared 



divergence in the soft endpoint region 
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(13) 



The typical values mg = (500 ± 200) MeV, and A = Aqcd='^50 MeV will be taken in our 
numerical calculations. Use of this gluon propagator is supported by the lattice J2J] and the 
field theoretical studies 2^, which have shown that the gluon propagator is not divergent as 
last as T^ . 

After giving the treatments for the infrared divergence arising from the gluon propagator, 
we now turn to show how to deal with a physical-region singularity of the on mass-shell quark 
propagators. It can be easily checked that each Feynman diagram contributing to a given matrix 
element is entirely real unless some denominators vanish with a physical-region singularity, so 
that the ie prescription for treating the poles becomes relevant. In other words, a Feynman 
diagram will yield an imaginary part for decay amplitudes when the virtual particles in the 
diagram become on mass-shell, thus the diagram may be considered as a genuine physical 
process. The Cutkosky rules J27| give a compact expression for the discontinuity across the cut 
arising from a physical-region singularity. When applying the Cutkosky rules to deal with a 
physical-region singularity of quark propagators, the following formula holds 

-k2-ks)% (14) 
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Xvb- 
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-ml_ 



-i7r6[{pb - k2 - kaf - ml], (15) 



where P denotes the principle-value prescription. The role of the 6 function is to put the 
particles corresponding to the intermediate state on their positive energy mass shell, so that in 
the physical region, the individual diagrams satisfy the unitarity condition. Eqs. (fT^ and p5p 
will be applied to the quark propagators Dd and Dt in Eqs. (| ^ -(| ^ . respectively. It is then 
seen that the big imaginary parts arise from the virtual quarks (d, b) across their mass shells as 
physical-region singularities. In fact, the above imaginary parts are among the main sources of 



strong phases for the B^ -^ K^K ,K*^K'^,K*^K* decays as discussed in the perturbative 



QCD approach 
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IV. DECAY AMPLITUDES AND LARGE STRONG PHASES 

Using the relevant input parameters hsted in Appendix A, we can calculate the tree and 
the penguin annihilation amplitudes for each decay mode and their corresponding numerical 
results for the quantities z and 6, which are presented in Table Q] With these considerations, 
the CP-averaged branching ratios for these decay modes are given in Table ITTl In this table, we 
present our "default results" along with detailed error estimates corresponding to the different 
theoretical uncertainties listed in Appendix A. The first error refers to the variation of the 
dynamical gluon mass, while the second one the uncertainty due to the CKM parameters 
A, A, p, and f/. Finally, the last error originate from the uncertainty due to the meson decay 
constants and the parameter fix- 

From the numerical results given in Tables. HI and ITT| we arrive at the following observations: 

• The strong phase associated with the tree amplitude At is about 6 ~ —45°, while the 
imaginary part of the penguin amplitude Ap is comparatively small; moreover, from the 
numerical results of the magnitude z, we can see that the penguin annihilation amplitudes 
are comparable to the tree ones, so that a large interference effect between the tree and the 
penguin annihilation amplitudes occurs. Combining these two ingredients, it is expected 

TABLE L The tree annihilation ampHtudes At, the penguin annihilation amplitudes Ap, magnitudes 
of z and the strong phase 6 for B^{B^) — > K^ K~ , K K^ , K*~^ K*~ decays. Here we give only the 
case when rUg = 500 MeV. 



Decay mode 


At 


Ap 


Z 


5 


SO(SO) - 


^ K+K- 


0.0273 - 0.032H 


-0.0114 + 0.0045i 


0.65 


-153° 


B^{B^) - 


-> K-K*+ 
-> K+K*- 


0.0467 - 0.0371Z 
0.0412 - 0.0464i 


0.0187 - 0.0005i 
0.0333 - 0.0092i 


0.70 
1.14 


37° 
33° 


B^{B^) - 


^ K*+K*~ 


0.0634 - 0.0543'i 


-0.0168 + 0.0047i 


0.47 


-155° 



TABLE II: The CP-averaged branching ratios (in units of 10"^) of B^{B^) -^ K+K~, K*^K^, 
K*^K*~ decays. The theoretical errors shown from left to right correspond to the uncertainties 
referred to as "gluon mass", "CKM parameters", and "decay constants and the parameter fix" as 
specified in the text. 



Decay mode 


Br 


Tt/T 


Exp. 


50(5°) - 


-* K+K- 


r, m +3.21 +0.83 +1.01 
■^•U^ -0.81 -0.54 -0.44 


- 


(4 ±15 ±8) X 10-s [21J 


B\B0) - 


-^ K-K*+ 
-^ K+K*- 


cr 7n +2.08 +2.26 + 1.11 
"J- '^-3.01 -1.71 -0.92 

A no +2.29 +1.63 +0.94 
^.zo_2.05 -1.21 -0.81 


— 


— 


BO(^O) - 


-* K*+K*- 


a OQ +1.78 +2.87 +2.04 
"•"^-1.50-2.19-1.58 


0.01 


< 1.41 X 10-^ \M 



that there exist large CP violations in these decay modes, which will be shown below. 

For a given dynamical gluon mass, the CP-averaged branching ratios of these decay modes 
follow the pattern: 

i3(i?° -^ K*+K*~) > B{B^ -^ K*^K^) > B{B° -^ K+K-), (16) 

which is due to the larger vector-meson decay constant fx* > fx for one vector meson in 
the final state or the larger spin phase space available for two final-state vector mesons. 

For the B^ — »■ K^K~ decay, the present obtained result is consistent with the experi- 
mental data 2l|, and also in agreement with the one given in J2i|: B{B^ -^ K^K") = 
(0.013lQ;oo5to;oo5lQ;o^]') x 10"^, when considering the huge uncertainties caused in the 
treatment of infrared divergence /q — ^ Xa = (1 + f)Ae*'^)ln^. The present theo- 
retical errors mainly originate from the variation of dynamical gluon mass, while our 
theoretical result alleviates the dependence of input parameters when compared to the 
one given by Ref. J30|. 

Another significant error comes from the parameter fix, as the decay amplitudes include 
/i|^ factor when considering the twist-3 wave functions contributions. Moreover, the 



present central value is smaller than the one given by the pQCD [37[: Br{B -^ K^K~) = 
4.6 X 10^^, which has used a bigger parameter /ix- If we also choose /ix = 2.22 GeV just 
as the pQCD does, and a dynamical gluon mass with mg = 420 MeV, we can get the same 

10 



result as the pQCD prediction. Namely, if choosing smaller fix, to get the same result 
as the pQCD prediction, one has to choose a smaller mg. In general, with the parameter 
fiK fixed at the value 1.4, 1.8 and 2.2 GeV, we can get the same result as the pQCD by 
taking rUg to be nig = 340, 370, and 420 MeV, respectively. Anyway, our prediction is 
consistent with the pQCD result after taking into account the theoretical uncertainties. 

• For the S° — > K^K*^ decays, from the numerical results we can see that the main 
theoretical errors originate from dynamical gluon mass and CKM parameters. In addition, 
we predict that BlB'^ -^ K~^K*~) ^ B{B^ -^ K~K*~^). The central values are larger 
than the ones given in Ref. |3o|: B{B^ -^ K^K*^) = (OMAtlfol-S-oml) x 10-^ but 
they remain consistent with each other within the uncertainties. 

• For the B^ — > K*^K*~ decay, only an upper limit at 90% confidence level has been 
reported |2a| 

B{B^ -^ K*+K*-) < 1.41 X 10"^ (17) 

Obviously, the present numerical result is far below the experimental data. It is noted 
that the branching ratio of this decay channel is less sensitive to the dynamics gluon mass, 
and the theoretical errors mainly come from CKM parameters. It is also noted that 99% 
of the branching ratio comes from the longitudinal part, and the result is consistent with 
that given by the pQCD method 38| |. 

V. LARGE CP VIOLATION 

We now turn to discuss the CP asymmetries in the B^{B^) -^ K^K~ , K*^K'^, K*~^K*~ 
decays. As there are the big strong phases and large interference effects in these decay modes, 
large CP-violating asymmetries are expected. It is very reasonable to neglect the transverse 
contribution and consider only the longitudinal part in B^ -^ K*^K*^ decay, since the trans- 
verse polarization fraction provides only 1% contribution to the total branching ratio of this 
mode. Thus we can discuss the CP asymmetries in the K^K~ and K*^K*~ decay modes by 
the same manner. 

Using the relevant formulas provided in the previous section, we are able to calculate the 
CP violations in the 5° ^ K+K' , K*^K'^, K*+K*- decays. Firstly, in Table ITTTl we present 
our predictions for the rephase-invariant CP-violating observables and their theoretical errors 

11 



in these decays. Secondly, taking the CKM angle a as a free parameter, the dependence of 
CP- violating parameters on the angle a is shown in Figs. El and |3] 

TABLE III: The CP asymmetries for 5°(fi°) -^ K+K' ,K*^K^ ,K*+K*- decays with the same 
error resources as in Table Hll 



mode 



Ckk 



ACkk 



Skk 



A.Skk 



K+K- 


n on +0.00 +0.04 +0.01 
"-"■"jy -0.04 -0.04 -0.01 


0.00 


n OR +0.09+ 0.04 +0.02 
U-OD -0.05 -0.08 -0.05 


K^K*^ 


n KQ +0.11 +0.05 +0.01 

"•^y -0.00 -0.07 -0.01 


n ni +0.00+0.01+0.01 

U-UJ^ -0.04 -0.01 -0.01 


n 7/1 +0.29+0.02+0.06 
*J-'* -0.12 -0.08 -0.04 


K*+K*~ 


n on +0.02 +0.05 +0.00 
"•<*'-' -0.01 -0.05 -0.00 


0.00 


n 70+0.01+0.09+0.00 
U-'O -0.01 -0.10 -0.00 



0.00 



n n7 +0.08 +0.06 +0.06 
-U.U( _o.08 -0.06 -0.06 



0.00 
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FIG. 3: The CP violation parameters Ckk and Skk for B^{B^) -^ K^K ,K*^K* decays as 
functions of the weak phase a (in degree). The dash-dotted, solid, and dashed lines correspond to 



rUg = 300 MeV, 500 MeV, and 700 MeV, respectively 



From Table UTTl and Figs. 3-4, we come to the following observations: 

(i) For the B^ —* K^K~ decay, due to the large strong phase 5 and large magnitude 
z ~ 0.65, the direct and mixing-induced CP-violating parameters Ckk and Skk are 
found to be quite large. The direct CP asymmetry is consistent with the one given 
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FIG. 4: The same as Fig.Elbut for S°(S°) -^ K^K*^ decays. 



by pQCD [37[. It is also seen that the CP asymmetry parameters Ckk and Skk are 
not sensitive to the choice of the dynamical gluon mass, and the main theoretical errors 
originate from CKM parameters. 

(ii) For the B^ -^ K*^K*~ decay, as it is dominated by the longitudinal part, the CP 
asymmetries in this decay mode have the same manner as the one in the B^ -^ K'^K~ 
decay. Thus we arrive at the similar conclusions as the ones for the B^ —>■ K^K~ decay. 

(iii) In contrast to the B^ -^ K^K~ , K*^K*" decays, the strong phases 6 in the B^ — >■ 
K^K*^ decay modes have opposite signs, so that the sign of the CP asymmetry parameter 
Ckk is also opposite to the ones in the K~^ K~ , K*~^ K*" decay modes. In addition, the 
rephase-invariant parameters Ckk and Skk in the i?° -^ K^K*^ decay modes also 
characterize large CP violation as the parameters ACkk and ASkk are small. Note that 
only the mixing- induced CP violation Skk is sensitive to the dynamical gluon mass. 

(iv) It is seen from figures |21 and E] that all CP- violating parameters have a strong depen- 
dence on the weak angle a. So these observables may be used to determine the range of 
the angle a in future more precise experiments. 
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VI. CONCLUSIONS 

In summary, we have calculated the strong phases, the CP-averaged branching ratios, and 
the CP asymmetries for the pure annihilation decays B^ -^ K~^K~ ,K^K*^, and K*^K*~ 
within the standard model. It has been shown that the nonfactorizable contributions dominate 
these decays, and the contributions of the penguin diagrams are comparable to that of the W- 
exchanged diagrams; the CP-averaged branching ratios of these decay modes are at the order 
of 10~^ ~ 10^^, and these small branching ratios predicted in the SM make them sensitive 
to any new physics beyond the SM. Of particular, as there are big strong phases and large 
interference effects between the tree and the penguin annihilation amplitudes, the CP-violating 
parameters Ckk and Skk have been predicted to be large in these decay modes. It has been 
seen that the CP-violating parameters have a strong dependence on the weak phase a, but they 
are not sensitive to the dynamical gluon mass except the mixing-induced CP violation in the 
B^ -^ K*~^K~ ,K^K*~ decays. The resulting branching ratio Br{B^ — ^ K'^K~) is consistent 
with the current experimental data. It is then expected that the predicted CP asymmetries 
should be reasonable. 

In this paper, we have adopted the Cornwall prescription for the gluon propagator with a 
dynamical mass to avoid the endpoint infrared divergence. Note that when the intrinsic mass 
is appropriately introduced, it may not spoil the gauge symmetry as shown recently in the 
symmetry-preserving loop regularization J39|. Meanwhile, we have also applied the Cutkosky 
rules to deal with the physical-region singularity of the on-mass-shell quark propagators. As a 
consequence, it produces the big imaginary parts which are very important for understanding 
large CP violations. The combination of the two treatments for the endpoint infrared divergence 
of gluon propagator and the physical-region singularity of the quark propagators enables us to 
obtain, by using the QCD factorization approach [l7|, reasonable results which are consistent 



with the ones 



37 



38j obtained by using the perturbative QCD approach[18[. However, different 
from the treatment of perturbative QCD approach, where the Sudakov factors have been used 
to avoid the endpoint divergence. As a consequence, it was shown that the pQCD predictions 
are insensitive to the choice of Sudakov factors and to the dependence of the impact parameter 
b U. 

It is noted that our present predictions for branching ratios depend on the dynamics gluon 
mass which plays the role of IR cut-off, such a dependence should in general be matched from 
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the nonperturbative effects in the transition wave functions. However this is not the case 
in general for direct CP violation. With the experimental and theoretical improvements, this 
quantity could be fitted from a well measured pure annihilation decay mode, and then expanded 
to other decays. In order to check the validity of the gluon-mass regulator method adopted 
to deal with the endpoint divergence, we plan to extend this method to the vertex corrections 
and hard spectator interactions for other charmless i?-meson decays. We expect that these 
corrections are independence of the dynamical gluon mass, which is being under investigation. 
Anyhow, the treatment presented in this paper may enhance its predictive power for analyzing 
the charmless non-leptonic 5-meson decays. 
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Appendix A: Input parameters 

The theoretical predictions in our calculations depend on many input parameters, such as 
the Wilson coefficient functions, the CKM matrix elements, the hadronic parameters, and so 
on. We present all the relevant input parameters as follows. 

The next-to-leading order Wilson coefficient functions (at /i = 772^/2) in the NDR scheme 
for 75 J22| have been used with the following numerical values 

Ci = 1.130, C2 = -0.274, C3 = 0.021, C4 = -0.048, 

C5 = 0.010, Ce = -0.061, C-j = -0.005/128, Cg = 0.086/128, (18) 

C9 = -1.419/128, Cio = 0.383/128. 

For the B meson wave function, we have taken the following results 



<^fip) = NBp\l-p)'exp 



1 f pniB X 

'2 V ^B ; 



(19) 



with ujb = 0.25 GeV, and Nb being a normalization constant. For the light cone wave functions 
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of light mesons, we use the asymptotic form as given in Refs. |41 , 







42 



$^(u) = $11 (n) = $^(n) = 6uu, twist - 2 LCDAs, 
0o-(n) = gl {u) = h\f\u) = Quu, 

h\^\u) = 3 - 12uu, 0fc(u) = 1, twist -3 LCDAs, (20) 

with u = 1 — u. 

For the other parameters used in our calculations, we list as follows J36| 

Ms, = 5.28GeV, mb = 4.66GeV, Mk'± = 0.89GeV, r^o = 1.536ps, 

Jb, = 200 ± aOMeV, Jk = 160MeV, fx* = 218 ± 4MeV, f^, = 175 ± 25MeV, (21) 

Vud = l- AV2, Vub = AX'ip - 17]), Vtd = A\%1 -p- 17]), Vtb = 1. 

The Wolfenstein parameters of the CKM matrix elements are taken as J36|: A = 0.8533 ± 
0.0512, A = 0.2200±0.0026, p = 0.20±0.09, r/ = 0.33±0.05, with p = p(l- f ), r/ = r/(l-^). 
The coefficient of the twist-3 distribution amplitude of the pseudoscalar K meson is chosen as 



;x^ = /i^ = 1.5 ± 0.2 GeV lli 



2^ 



Appendix B: The decay amplitudes 

To evaluate the hadronic matrix elements, we may adopt the QCD factorization formalism. 
For the annihilation process B -^ M1M2, the matrix element can be written as J30| 

(MiMalO.I^O) = fB^B * /a/,$a/, * Im.^m, * T„ (22) 

where Oi is the effective operator appearing in the effective weak Hamiltonian, the * products 
imply integrations over the light-cone momentum fractions of the constituent quarks inside 
the relevant mesons. Tj is the hard-scattering kernel that can be computed perturbatively 
with the QCD factorization approach. $jv/ and /« are the leading-twist light-cone distribution 
amplitudes and the decay constants, respectively. 



Using the formula ()22|) and the meson wave functions given in J40|, |4l|, |42, |43[, we can 
evaluate the decay amplitudes of the PF-exchanged diagrams in Fig. ^ (only the Oi operator 
has contribution) and the penguin annihilation diagrams in Fig. |21 (O4, Oq, Og, and Oio have 
contributions). 
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For final states containing two pseudoscalar K meson, the amplitudes are found to be 



wliere 



AT{Bd -^ K+K-) = Ci X Ai, 



2C, + ^jxA,+ (2C, + ^]xA2, 



A, = ^f^fl^as{ii)% f di f dx ['dy<^nO 
V2 ^^c JO Jo Jo 

( f /l/f 4 /l#4 \ 



DA'' 



D.k^ 



+ 



Mil 



[x + y- i)(l)k{x)(t)k{y) + {y~i- x)(l)k{x 



6 



(y-^-a;)0fe(y) 



0;(x) 



{x + y-0 



<P'Ax)4>'Ay)\ M] 



Kiy) 



6 ' " '' 36 

+ Ux + y + ^ - 2^B)Mx)My) + (y + ^ - x)Mx) g 



(23) 



(24) 



A, 



X 



+ 



^/i./ivra.(/i)^^ d^J^ dxj^ dy$f(0 

{y-i)T7T-. + i^-^B)^MK{x)^K{y) 



DA'' 



D,k' 



i4 



{x + y- Ci<Pk{x)(t)k{y) -{y-i- x)4)k{x 



^ €iy) 

6 



+ yy-i- x)My)^ [x + y-O- 



6 ' " "' 36 

{x + y + ^- 2^B)Mx)My) -{y + ^- x)^^) 



Dak' 

<P'Ay) 



6 



+ {y + i - x)(t)k{y) 



<p',{x) 



+ {x + y + i-2) 



<P',{xWM\ M 



'B 



DA'' 



(25) 



6 36 / lyb^ 

wliere ^b = {Mb — nib) /Mb with Mb being the mass of Bd meson, (j)„{x) = ^^^ , $'s and 0's 
are the leading twist (twist-two) and twist-three light-cone distribution amplitudes of mesons, 
respectively. We set the scale /i to be the averaged virtuality of the time-like gluon, i.e., 
/i = mb/2. k' and Db^d arise from the propagators of the virtual gluon, the bottom quark 6, 
and the down quark rf, respectively. 

For Bd — > K^K*^ decays, we need only consider the longitudinal wave function of the 
vector K* meson, due to the conservation of the angular momentum. The corresponding decay 
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amplitudes of Bd -^ K^K* are given with the same form as Eq. (f^ . but with the amphtudes 
of Ai and A2 replaced by 



A, = ^f^f^nasifi)^ f d^ f dx f dy^fiO 

V ^ ^^C JO JO Jo 



X 






+ fk-'f^K 



rriK* 



Mb 
+ {x + y-0 



(-x + y-0 



hi'\yWAx) 
12 



[x + y-0 



ht\ymx) 



^ + 2{x-y + 0<P{x)h\\y)^ ^ 



I)^fc2 



Ms 



-x + y-0 



h'l'\y)4>'Ax) 
12 



(x + y + e - 2e 



/^?H2/)0(^) 



B 



-i + (x - y - e)0(a;)/i}, '(y) ^ "^ 



/^/,A;2 



(26) 



A, 



^-jBfKTrasif^)^^ I di j dx I dy<l>f{^) 



V2^ 



X 



/. 



K* 



NcJo 
^ Mi , 

-^ + ^)z^ + ^^ + ^-^ 






+ fK*(^K 



:(-), 



Mb 

- {x + y-0 

. rriK* 
+ ^^'''^ M^ 

+ ix + y + ^-2) 



<ih\^\y) 



(-x + y-0 



hl'\y)4>'Ax) 

12 



$i^(x)$i^(i/) 



+ (x + 2/ - 



ht\ymx) 



^ + 2{x-y + i)<p{x)hfMy) ^ ^ 



/^^A;2 



(-X + y - 0^ T7^ + (x + 2/ + e - 2^^)^ ^ 



12 






/^/,A;2 



(27) 



where /in (y) = — ^ — . With the change for the signs of the second and the third terms in the 
twist-3 amplitudes in Eqs. ()2fi|l and (j2Zj) and the exchange of the variable x and y of the wave 
functions, we can get the corresponding decay amplitudes of B^ -^ K~K*^ mode. 

Finally, for the B — * K*^K*~ decay, in the rest frame of B^ system, we have Ai = A2 = 
A (where Ai and A2 denote the helicities of the K*^ and K*~ mesons, respectively) through 
the helicity conservation since the B^ meson has helicity zero. So, there are generally three 
decay amplitudes, Hq, H^, and H_, representing A = 0, 1, and —1, respectively. Considering 
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only the leading twist contributions, we can get the longitudinal amplitudes Aot and Aop of 
the decay B^ — > K*^K*~ from Eqs. (j23) and (j^^ by keeping only the $i^(x)$j^(|/) term and 
replacing the decay constant fx by fx*- The total longitudinal amplitude is then given as 

^0 = KbKTdAT - Vtiy:^Aop. (28) 

As for the transverse amplitude if±, we have 

+ {-f{x)g^l^{y) + g^l\x)gt\y)/A - f{x)g}_^\y)/8 + (7>)(x)^i:)(y)/32) ^^ 



^/./i^^-.(.)| /^ ^^ /' ^^ /^ ^^^^^^) (2^^ + ^ 



X (^(-^f (x)/(y) + ^f (x)^?? )(y)/4 - gT\x)fiy)/8 + ^(_«)(x)^>)(y)/32)-0 

+ (-^?^(x)/(y) -#(x)^f)(y)/4 + ^>)(x)/(y)/8 + (7i")(x)^>)(y)/32)-^ (29) 

where the function f{x) = j^ du (0||(m) — gl {u)\ and g}_ (y) = ^\ ■ Changing the signs 
of the g_l^ gj^ and g}^ f terms in Eq. pH)) . we can get the other amplitude H^. 
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